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Abstract 

A local void in the globally Friedmann-Robertson- Walker cosmological model with the critical 
density (Qo = 1) is studied. The inhomogeneity is described using a Lemaitre-Tolman-Bondi 
solution for a spherically symmetric distribution of matter. The scale of the central underdense 
region is ~ 150 Mpc. We investigate the effects this has on the cosmological time scale, the 
measurement of the Hubble constant and the redshift -luminosity distance for moderately and very 
distant objects (z ~ 0.1 and more). The results indicate that if we happened to live in such a void, 
but insisted on interpreting cosmological observations through the FRW model, we could go wrong 
in a few instances. For example, the Hubble constant measurement could give results depending 
on the separation of the source and the observer, the quasars could be younger than we think and 
also less distant (less energetic). 
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1 Introduction 



It seems, particularly after the introduction of the inflationary paradigm ||J], that the isotropic and 
homogeneous Priedmann-Robertson- Walker (FRW) cosmological models are best suited for the 
description of the global structure and the evolution of the universe. However, a similar statement 
is not necessarily true when cosmologically moderate scales are thought of. 

There exists direct observational evidence in favour of the isotropy of the observed universe, 
namely, the COBE data confirming a high degree of isotropy of the cosmic microwave background 
radiation (CMBR) [||, gj. However, there is no observationally based reason supporting the as- 
sumption of homogeneity. On the contrary, there seems to exist observational evidence in favour 
of larger and larger structures P] . 

Recent work on modelling voids in the expanding universe ||] has shown that it is possible to 
construct an asymptotically FRW universe containing an expanding spherical void of the Lemaitre- 
Tolman-Bondi type i. Moreover, the void can be a region of under-density rather than vacuum. 

We think that there exists sufficient observational evidence (briefly discussed later in this 
paper) to support a conjecture that we may live in a relatively large underdense region embedded 
in a globally FRW universe. Exploring physical properties of such a model is the aim of the present 
paper. 

In the following section, we briefly discuss the LTB model. Section || consists of a brief 
discussion of the observational background and the simple toy model of a local void presented 
here. The closing section contains a description of our results of numerical calculations as well as 
conclusions. 

Throughout this paper we use units in which G = c = 1, unless stated otherwise. 



2 The Model 



First, for the sake of notational clarity, let us recall the FRW line element: 

dr 2 



ds 2 = dt 2 - a 2 (t) 



1 — kr 2 



+ r 2 dn 2 



(1) 



with dQ 2 = d9 2 + sm 2 6d(b 2 . 



1 A cosmological solution spherically symmetric about one point was first proposed by Lemaitre However, it 
is usually called the Tolman-Bondi solution Jt| [s| . 
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Now, let us consider a Lemaitre-Tolman-Bondi |?], |8| model for a spherically symmetric in- 
homogeneous universe filled with dust. The line element in comoving coordinates can be written 
as: 

ds 2 = dt 2 - R' 2 (t, r)f- 2 dr 2 - R 2 (t, r)dCl 2 , (2) 
where / is an arbitrary function of r only, and the field equations demand that R(t,r) satisfies: 

2RR 2 + 212(1 - f 2 ) = F(r), (3) 

with i* 1 being an arbitrary function of class C 2 , R = dR/dt and R' = dR/dr. We have three distinct 
solutions depending on whether f 2 < 1, = 1, > 1 and they correspond to elliptic (closed), parabolic 
(flat) and hyperbolic (open) cases, respectively. 
The proper density can be expressed as: 

F' 



H WttR'R 2 ' 

Whatever the curvature, the total mass within comoving radius r is: 



(4) 



4 

so that 



M(r) = - / drf^F' = 4tt I drpf^R'R 2 , (5) 



o J o 



M'(r) = ^=47rpf- 1 R'R 2 . 
dr 

Also for p > everywhere we have F' > and R' > so that in the non-singular part of the model 
R > except for r = and F(r) is non-negative and monotonically increasing for r > 0. This 
could be used to define the new radial coordinate f 3 = M{r) and find the parametric solutions for 
the rate of expansion. 

In the flat (parabolic) case f 2 = 1, we have 

R=l(9Ff 3 (t + (3f 3 , (6) 

with (3(r) being an arbitrary function of class C 2 for all r. After the change of coordinates R(t, f) = 
f{t + /3(f)) 2 / 3 , the metric becomes: 

ds ^ = dt 2 - (t + f3f 3 (Y 2 dr 2 + r 2 dQ 2 ) , (7) 

where 
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and from (Q) the density is given by 

P = (9) 

Clearly, we have that (t — > oo) the model tends to the flat Einstein-de Sitter case. 

For the closed and open cases the parametric solutions for the rate of expansion can be written 
as §: 

i? = ^(l-/ 2 )~ 1 [l-cos( ? ;)] 5 f<l, (10a) 
t + (3=±F(l-f 2 y V2 [v-sm(v)], / 2 <1, (10b) 



and 

R 



\F(f 2 -l) 1 [cosln» - 1] , / 2 >1, (11a) 
t + (3=^F(f 2 -iy V2 [ S mh(v)-v], / 2 >1, (lib) 

with P(r) being again a function of integration of class C 2 and v the parameter. 

The flat case (f 2 = 1) has been rather extensively studied elsewhere ||10|| . The model depends 
on one arbitrary function f3(r) and could be specified by assuming the density on some space-like 
hypersurface, say t = to- 

The cases of interest to us, ( |i~0| ) and (p"l"D, correspond to closed and open models, respectively. 

Before we proceed (in the next section) to discuss the observational grounds for modelling a 
local void, we need to amplify the discussion of the LTB model by introducing basic features of the 
propagation of light in our model. The high degree of isotropy of the microwave background forces 
us to the conclusion that we must be located very close to the spatial centre of the void. In our 
discussion, for the sake of simplicity, we place an observer at the centre (tob = to-, fob = 0). 

The luminosity distance between an observer at the origin of our coordinate system (to, 0) and 
the source at (t e ,r e ,9 e ,(j) e ) is ||: 

dL= \Mf) = R ( t e,re)[l + z(t e ,r e )} 2 , (12) 

where £ is the absolute luminosity of the source (the energy emitted per unit time in the source's 
rest frame), J- is the measured flux (the energy per unit time per unit area as measured by the 
observer) and z(t e ,r e ) is the redshift (blueshift) for a light ray emitted at (t e ,r e ) and observed at 
(i„,0). 
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The light ray travelling inwards to the centre satisfies: 



ds 2 = dt z - R' 2 (t, r)f-' z dr z = 0, d6 



0, 



and thus 



^- = -R'(t,r)/f(r). 



Without getting into a detailed discussion, which can be found in 
the equation of the light ray travelling along the light cone is: 

t = T(r), 

using (|l3|) we get the equation of a ray along the path: 



(13) 



let us state that if 



dT{r) 



dr 



R' 

-y[ T ( r ), r ] 



(14) 



(15) 



where 



R'[T(r),r] 



d 2 R 



dtdr 



r,T(r) 



am 

~~di~ 



r,T(r) 



The equation for the redshift considered as a function of r along the light cone is: 

dz 

- = (l + z )R'[T(r),r], 

and the shift z\ for a light ray travelling from (ii,ri) to (to,0) is: 

M'(r) 



log(l + zA 



log(l - Ol) 



dr 



o r(l-oi)' 



(16) 



(17) 



where 



ai (r) = R[T(r),r], 



and, in obtaining the second equation, we used (^) and (||). Thus we have two contributions to the 
redshift. The cosmological redshift due to expansion, described by the first term with a\ = R, and 
the gravitational shift due to the difference between the potential energy per unit mass at the source 
and at the observer. Obviously, in the homogeneous case (M'(r) = 0) there is no gravitational shift. 



3 The modelling of the local void 

If we restrict ourselves to spatial scales that have been well probed observationally, i.e. up to a few 
hundred Mpc, the most striking feature of the luminous matter distribution is the existence of large 



4 



voids surrounded by sheet-like structures containing galaxies (e.g. The surveys [11], [12] give a 



typical size of the voids of the order 50-60 Mpc. There has also been some evidence [13] -with 
less certainty- for the existence of larger underdense regions with characteristic sizes of about 130 
h^ 1 Mpc. Also, dynamical estimates of the FRW density parameter give very different results on 
different scales. The observations of galactic halos on scales less than about 10 to 30 Mpc typically 
give (see e.g. ]14|| ) O10-30 — 0.2 ± 0.1. On the other hand, smoothing the observations over larger 
scales (> 20 Mpc, say ~ 100 Mpc ) indicates (e.g. [fl2|| ) the existence of a less clustered component 
with a contribution exceeding 0.2, and perhaps as high as fi^ioo — 0-8 ± 0.2. 

At the same time, the large scale galaxy surveys (some of the recent literature is given in 



[15]) firmly indicate a considerable excess in the number-magnitude counts for faint galaxies rel- 
ative to predictions of homogeneous, "no-evolution" models. This excess could be the result of 
a non-standard galactic evolution or could be caused by rather exotic FRW cosmology (i.e. the 
deceleration parameter qo <C 0.5 or a non-zero cosmological constant A). However, it can also be 
treated as an observational indication of a very large (on the scale of the redshift z ~ 0.5) void. 
A model of a local void with a density distribution based on the faint galaxies number counts is 
presently being studied |i~6f |. 

In the model presented here, we confine ourselves to the simple density distributions. We 
study two cases: a void with the central density equal to that of an FRW model with the density 
parameter Qq = 0.2, asymptotically approaching the FRW model with £Iq = 1, and a very similar 
void "distorted" on an intermediate scale by a peak (O ~ 1) in the density distribution. The two 
distributions are: 



and 



^ 1 r \ 2 exp(r/L) 



LJ \exp{r/L) — 1] 



(18a) 



r\ exp(r/L) 



+ ( j)V('/'> a . (.181")) 



,LJ [exp(r/L) — l] 2 

In the numerical calculations presented in the next section we used the values £l m in = 0.2, Vl max = 1 
and L = I = 30Mpc. This assures that the void converges satisfactorily fast to the outside critical 
FRW universe (O ~ 0.86 for r ~ 150Mpc and ~ 0.95 for r ~ 200Mpc) and that the intermediate 
peak is observationally acceptable (f2 ~ 0.64 for r ~ 33Mpc ). 
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4 The results and discussion 



In general, an LTB model depends on three arbitrary functions, see section ^, F(r),(3(r) and f(r). 
Since F(r) can be interpreted as twice the effective gravitational mass within comoving radius 
r ([§]), then, in accordance with the discussion following (|5|), assuming its form is equivalent to 
a coordinate choice. In our calculations we used F(r) = 4r 3 . The second function, /3(r), sets 
the initial singularity ("big bang") hypersurface of the model. Since we want the outside region 
in our toy model to be fully equivalent to the critical FRW universe, we set (3{r) = 0, thereby 
assuming a universally simultaneous big bang. We also set the time coordinate of constant time 
hypersurface "now" so that it is equal to the age of the universe to in the FRW model with Qq = 1. 
In doing so, we give up a very important feature of an LTB model: an extra (with respect to 
FRW) degree of freedom that would allow the age of the universe to be different from FRW or even 
position dependent. The third ("curvature") function, /(r), is an unknown to be solved for in our 
calculations. From the work done in ||, we conclude that the LTB case to be used in modelling an 
underdense comoving void in an FRW universe is the hyperbolic (/ 2 > 1) one. 



In a manner similar to that employed in [1C], we assume that since all cosmological observa- 
tions are necessarily done by detecting some form of electromagnetic radiation, the solution should 



progress along the light cone. The final set of equations we solve consists of the equations (15), 
( |i~(f ) and the equation describing the density distribution (||) through either of the relations (18) 
taken along the light cone, e.g.: 

r > = ifrmnWr] = n " (r) - (19> 

Since T(f) (the time of emission t e of a light ray observed at r = at to) is now given by ( |1 lb|) , 
the functions to be solved for are f(r), z[r) and v(r), where the parameter v becomes the function 
of position. The initial conditions for the integration have to be set at r ^ 0, since the analytic 
expressions (11) are singular at r = 0, where f 2 = 1 (we have a flat Qo = 0.2 FRW universe 
there). We assume that for the initial radius rj <C 1 (we use dimensionless radius and time in the 
calculations) corresponding z% and tj are given by their standard FRW values. Then z(ri),v(ri) 
and can be obtained from ( |llb[ ), ( |l5| ) and (|i~6|). 

Once the equations have been numerically integrated we use ( llbj ) to obtain t e for a given r e . 
The luminosity distance di corresponding to this event is obtained with the use of ( |I2"D and z{t) 
(useful in studying the cosmological time scale) is given by the parametric relation [T(r), z(r)]. 



G 



FRW 0.2 

00 0.2 0.4 0.6 r 0.8 1 1.2 

Figure 1: The redshift z as a function of comoving radius r. The FRW results for ft = 1 and fl = 0.2 
are denoted by "FRW 1" and "FRW 0.2" , respectively. The results for the LTB modelled voids are denoted 
"voids" . 

The results of our numerical calculations are as follows. Figure || depicts z(r), where r is 
the dimensionless comoving radius used in the calculations. The coordinate distance has no direct 
physical relevance, but our units here are such that r = 1 corresponds to 2997. 95/i -1 Mpc, where 
h is the usual coefficient in the observationally determined value of the Hubble constant: Hq = 
WOh km s^Mpc -1 . 

The departure of z(r) from its FRW behaviour does not seem to be dramatic. In fact, on the 
distance scale used in Figure [j] it is hardly noticeable. More details, on a smaller spatial scale, can be 
seen in Figure ||[ The redshift z, after being influenced noticeably by rapid changes in the density p 
on smaller scales, asymptotically tends to a limit that could, in accordance with FRW interpretation, 
correspond to the universe with the density parameter in the range Vl € (0.2; 1). i The increase in 
z on intermediate scales is clearly induced by the additional gravitational shift caused by the mass 
distribution of the void. The large scale behaviour is controlled by our assumption of the equality 
*o = toFRW- 

2 This situation does not change considerably if we model the central void in a similar manner but with respect to 
the observable z as a variable H . 
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Figure 2: The redshift z as a function of comoving radius r. Smaller scales. Notation as in Fig. 1. 



However, one should not forget that the comoving distance is not an observable, whereas the 
luminosity distance is. In principle, provided we know its absolute luminosity C, we can establish 
the luminosity distance, defined by AttoIl 2 = C/J-, by measuring the energy flux J- of an observed 
object (for a discussion of usual caveats associated with so-called "standard candles" see e.g. [17|). 

Due to the lack of space we do not present the z(t) relation here. There is, again, some small 
scale (to — t -C 1) divergence from the FRW behaviour, but for early cosmological times (t <C 1) 
the relation tends to the critical FRW (O, = 1) one. This is in accordance with our assumption 
of a simultaneous big bang, (3(r) = 0, and with our setting the age of the universe to be equal to 
that of the critical FRW case. Objects with redshifts of order a few are younger than their FRW 
counterparts, but not significantly. 

Figures |3] and || present the redshift -luminosity distance relation on cosmologically large (Fig- 
ure |3|) and intermediate (Figure ||) scales. Now the departure from the FRW behaviour becomes 
apparent. This comes as no surprise if one recalls the formulae for di in both LTB dl2| ) and FRW: 



d L 2 = a 2 (t )r e 2 (l + z e 



(20) 



The important question is whether our results contradict the linearity of the Hubble relation z 
Hd, well established on small scales. 
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Figure 3: The rcdshift z as a function of the luminosity distance di (in Mpc). 
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Figure 4: The redshift z as a function of the luminosity distance di (in Mpc). Smaller scales. 
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Due to our choice of the cosmological time scale and FRW embedding, the asymptotic (z — > oo) 
behaviour of z(d£) is that of the critical (O = 1) FRW case. On intermediate scales, however, 
objects of comparable redshifts are located at smaller luminosity distances. The ratio diTB/dFRW 
is « 0.7 (« 0.5) for z « 1, « 0.8 (« 0.4) for z « 2.5 and « 0.9 (« 0.36) for z « 4.5, where the most 
distant quasars are observed. (Values in parentheses correspond to the FRW O = 0.2 case.) Since 
the absolute luminosity C of the source scales as the square of dt this reduces the energy output of 
QSOs (up to an order of magnitude). Also, the angular diameter distance dj± = D/5 = (1 + z)~ 2 dL 
now gives a smaller proper distance D across the source for the same observed angular diameter 
5. This helps resolve problems with seemingly acausal signals (correlations in luminosity bursts) 
observed across some quasars. 

At the same time, inspection of Figure |4] shows that on small scales a very nearly linear (in fact, 
observationally indistinguishable from linear) "Hubble diagram" is obtained. However, a different 
value for the Hubble parameter (constant) is inferred (position, or rather d^, dependent on larger 
scales), if we insist on interpreting the results of cosmological observations through an FRW model. 

To explore this possibility let us recall that in FRW cosmology the exact result for the Hubble 



relation {z versus d£) in the matter dominated universe is [17] 



H d L = q - 2 [zq + (q - 1) [^2zq + 1 - lj j , (21) 

where go = —d(to) / a(to) Ho 2 is the deceleration parameter. 

Let us assume that we live in a local LTB void and the z vs. d^ relation differs from the 
FRW one as described in this paper, but we are biased by our theoretical prejudice and interpret 
cosmological observations through the FRW model. 

On cosmologically small distances we measure the same value of Hq independently of the 
model (we call this value "the local measurement"). This stems from the fact that, due to our 
assumptions, very close to the centre (r<l) the model is well approximated by the FRW universe 
with f2 = 0.2. Obviously, if the universe were LTB rather than FRW, then the Hubble parameter 
based on the observed (LTB) values of z and c^, but inferred through an FRW relation fl2l|), would 
be position (redshift) dependent as shown in Figure ||. 

The values of Hq reported to date span the range 40 to 100 km s~ 1 Mpc _1 (with standard 
errors quoted frequently as 10 km s _1 Mpc _1 or less!). Inhomogeneities similar to the LTB void 
presented here might provide an explanation for this. 
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Figure 5: The "observed" Hubble constant Hq (in units of the local measurement) as a function of the 
redshift z. The results interpreted through FRW ft = 0.2 and ft = 1 models, respectively, are denoted "int. 
as FRW 0.2" and "int. as FRW 1". 

The LTB void, such as the one presented here, decreases its density contrast (the depth of the 
void with respect to the FRW background) when evolved back in time 10|. At early times it is 
almost homogenized (at t/to ~ 10~ 5 we have \pltb(t)i ' Pfrw — 1| < 10 -6 ). This corresponds to a 
universe which at the beginning is very similar to the FRW one, but different at late stages. 

In this manner, while retaining all accomplishments of the FRW cosmology in dealing with 
epochs preceding the matter dominated era, we can gain new freedom in modelling the more recent 
universe. We can solve the age of the universe problem (by assuming (3{r) = const ^ 0), provide 



the excess power observed on scales of 5 — 10,000 km s _1 in modelling structure formation (see [10]), 
alleviate a few old problems associated with quasars (their age, luminosity and size) and provide 
an explanation for the wide range of reported values of the Hubble constant. 
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